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ABSTRACT: 


it  study  of  the  statistical  theory  of  light  propagation  in 
a turbulent  medium.  A parabolic  equation  which  was  examined 
was  only  valid  when  back  scattering  could  be  disregarded. 

The  theory  proposed  uses  the  small,  value  of  the  longitudinal 
acale  of  heterogeneities,  as  compared  with  other  longitu- 
dinal scales  of  the  problem.  Many  difficulties  encoun- 
tered in  the  range  of  strong  fluctuations  are  eliminated 
In  the  new  theory. 
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Introduction 

In  the  propagation  of  electromagnetic  waves,  in  a medium  with  random 
large-scale  (in  comparison  with  the  long  wave)  irregularities  due  to  the 
multiple  forward  scattering  effect,  fluctuations  of  the  wave  field  rapidly 
increase  with  distance.  Beginning  at  a certain  distance,  they  are  unsuit- 
able for  calculation  by  the  perturbation  theory  in  any  of  its  forms  (range 
of  strong  fluctuations)1.  This  effect  was  observed  experimentally  by 
Gracheva  and  Gurvich  in  experiments  on  the  propagation  of  light  in  a 
turbulent  atmosphere,  and  in  further  detailed  research  in  work  by  Gurvich, 

Kail Istratova,  Time  [3]  , Gracheva,  Gurvich,  Kellistratova  £4] , and  Mordu- 
kovioh  [5l. 

Recently,  a number  of  books  have  appeared  in  which  equations  have  been 
obtained  by  different  methods,  and  which  describe  the  strong  fluctuation  region 
of  a field  t6  - 1?J.  The  method  for  obtaining  these  equations,  used  In  works 
[10  - 1?]  , is  based  on  the  approximation  of  the  wave  propagation  process  in 
a nonuniform  medium  by  the  random  diffusion  process.  In  this  approximation, 

^The  range  of  fluctuations  of  intensity,  which  is  describod  by  the  first 
approximation  of  the  even  perturbations  method,  is  called  the  weak  fluc- 
tuations range.  The  basic  results  of  theoretical  and  experimental  research 
in  this  range  have  been  shown  in  some  detail  in  a book  by  Tatarskiy  [ll  (see 
also  below  Section  6 of  the  first  chapter). 
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closed  aquations,  suitable  for  the  strong  fluctuations  range, can  be  obtained 
for  all  moments  of  a field  and  equation  of  the  Binetein-Fokker  type  can  be 
obtained  for  the  characteristic  function  of  a field. 

In  this  work,  we  try  to  show  in  sequence  the  basio  results  obtained  in 
this  direction.  The  work  oonsists  of  two  chapters. 

The  first  chapter  is  devoted  to  a general  examination  of  the  propagation 
process  of  a light  wave  in  a turbulent  medium,  The  first  section  examines 
the  stochastic  equation,  desoribing  propagation  of  a wave  in  a random  medium 
and  several  aoourate  conclusions  of  this  equation  are  given.  The  seoond 
seotlon  examines  a model  in  which  one  can  disregard  the  longitudinal  radius 
of  the  dieleotrlc  constant  correlation  in  comparison  with  all  the  longitudinal 
scales  of  the  problem.  This  supposition,  equivalent  to  the  substitution  of 
an  actual  correlation  function  of  the  refractive  index  on  the  delta-function 
in  a longitudinal  direction,  allows  one  to  obtain  closed  integral  equations 
for  all  moments  of  the  wave  field,  Throughout  the  work,  apart  from  the 
supposition  above,  there  is  also  a supposition  on  the  Gauss  distribution  of 
probabilities  for  fluctuations  of  the  refractive  index.  Here,  one  can  reduce 
the  integral  equations  to  differential  ones  and  show  that  the  characteristic 
functional  of  a field  is  satisfied  by  an  equation  of  the  Einstein-Fokker  type, 
and  shows  that  the  propagation  of  a wave  is  a diffusion  process.  The  fourth 
section  shows  a method  of  successive  approximations  for  solving  the  stochastic 
equation  of  a wave  propagation,  in  which  the  diffusion  approximation  shown 
above  is  the  first  approximation.  Investigation  of  the  second  approximation 
allows  one  to  obtain  boundaries  for  the  applicability  of  the  diffusion  approxi- 
mation and  show  that  the  latter  can  also  be  used  in  the  strong  fluctuation 
range  of  a field.  In  the  fifth  section  there  is  an  examination  of  an  actual 
example,  devoted  to  calculating  the  mutual  coherence  function  in  a turbulent 
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medium  and  examining  the  comparison  of  results  of  calculation  with  experi- 
mental data.  The  final  section  of  the  first  chapter  investigates  the  ampli- 
tude-phase characteristics  of  a light  wave. 

In  a whole  series  of  works  describing  wave  propagation  in  a medium  with 
random  irregularities,  the  equation  for  ray  diffusion  Is  used  [18  - 2t)  (an 
approximation  of  geometrical  optics).  The  diffusion  equation  itself  (the 
Einstein- Fokker-EEF)  is  normally  written  on  the  basis  of  intuitive  considera- 
tion relying  on  the  analogy  with  well-known  problems,  leading  to  this  equation. 

The  dynamic  equation  of  the  problem  (in  this  case  — the  equation  of  rays)  is 
only  used  for  calculating  coefficients  in  the  EEF.  As  was  said  in  work  [22,  23), 
it  is  still  not  clear  in  what  conditions  one  can  justifiably  use  the  EEF  (a 
number  of  such  conditions  was  suggested  in  [23)  from  physical  considerations). 

At  the  same  time,  there  must  be  certain  conditions  superimposed  on  functions 
entering  the  dynamic  equations,  and  certain  limitations  on  parameters,  on 
which  the  solution  depends,  in  which  the  EEF  will  be  the  logical  result  of 
dynamic  equations.  These  conditions  are  examined  in  work  [24] , where  it  was 
shown  that  the  EEF  can  only  be  valid  in  narrow  angle  approximation.  Sta- 
tistical characteristics  of  amplitude  and  phase  of  a light  wave,  as  shown  in 
work  [25]  , are  determined  by  the  statistical  characteristic  of  rays. 

In  the  second  chapter  of  this  vwrk,  light  propagation  in  a turbulent 
medium  is  investigated  in  the  approximation  of  geometrical  optics.  The  first 
seotion  examines  the  problem  on  the  diffusion  of  rays  in  random  irregular 
media,  in  the  second  and  the  third  seot.ions  the  amplitude-phase  fluctuations 
of  a light  wave  are  examined, 
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I,  General  Examination 


1,  Initial  stoohastio  aquations  and  some  of  their  consequences 
The  propagation  of  mono  chromatic  light  in  a medium  with  large-scale 
irregularities,  when  depolarization  is  small  026],  can  be  quite  accurately 
described  by  the  scalar  wave  equation 

+ +V)]f  -0.  (1) 


Here  ft  is  linked  with  a component  of  the  electric  field  E by  the  ratio 


F.  *a  i exp(—  M),  -p  <«>,  • ws  * ^ — which  is  the  fluctuating  part 

C V / 


of  the  dielectric  constant.  If  one  disregards  large  angle  scatter,  then 
instead  of  (1),  one  can  use  the  parabolic  equation  for  function  u linked 


with  ^ by  the  ratio  ^ u «xp(ikx)i 


...  On 


2 /*  -■--  T AxM  -f  f)u(x,  p)  **  0, 

>-<>'•'1'  V-gr+'jF- 


(2) 


During  the  transition  from  (1)  to  (2),  the  term  d^/dx2  is  rejected.  The 
initial  condition  for  (2)  is 

p)  **  w#(p). 

Further,  we  shall  proceed  from  (2). 

Equation  (2)  with  the  boundary  condition  (3)  can  be  written  in  the  form 
of  an  integral-differential  equation 

X 

k 


u(x,  p)  •>  «n(p)cxp  cf  5 ^ (S,  p)J  4- 


(4) 


X x 

+ 2^,frf5eXp^2  p), 


which,  in  several  cases,  is  more  suitable. 

If  one  writes  the  solution  of  equation  (4)  in  the  form  of  an  interative 
aeries,  one  can  easily  be  persuaded  that  function  u(x(/i 0)  will  only  depend 
functionally  in  the  previous  values  along  x of  i(|  ,p)  from  the  interval 


^ 
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x,  It  follows  that  u(x,^)  doss  not  change  whan  varying  the  func- 
tion outside  of  this  Interval,  that  Is,  In  sections  £<0,  £ > x. 

Consequently,  the  variation  derivative  rfu/</)f  satisfies  the  oondltlon 


Jffo-lL  ~ o (*'  < 0,  .v'  > x), 

*«(*'.  P') 


(5) 


which  we  shall  call  the  causality  condition.  Furthermore,  a more  suitable 
&t(x,  p) 

value  for  us  is  ~ , “*  where  x = x1.  This  value  can  be  found  from 
W.  P ) 

equation  (2),  if  it  is  integrated  in  respeot  to  x within  limits  (0,  x), 

& 

and  afterwards  functioned  by  the  operator 


&*(*',  p') 


+ f P)1  + 

pO  V'  p’) 


+ A^(p  - P>(/,  f')  - o. 


Here  we  calculated  that  p)-  .v') 6(p  — p'),  and  the  lower  limit 


li(x\  p') 


for  Integrating  0 was  replaced  by  x* , since  according  to  (5)  «() 


when  *'>?,  . Assuming  that  x*  = x,  we  obtain  the  ratio  (10) 


.ij^L  p - p>(.v,  P). 


3s (x,  P') 


(6) 


_ , S u(x.  p) 

In  most  cases,  the  value  — --- 


when  .v  can  be  expressed  by  the 

CefV.  p'> 

funotion  of  the  Green  equation  (2),  linking  u(x, fi)  and  u(x'^0' ) when 


0 < x'  < xi 


II (X,  P)  “ (dp'O’ix,  Pi  .v\  ?')«(*',  p'). 


by  using  the  ratio  M > 


„ it  G(x.  p,  x',  p')«(x\  p'). 


*«(*'.  p') 

The  solution  of  equation  (2)  with  the  initial  condition  (3)  can,  by  using  the 
methods  suggested  by  Fradkin  [27,  28],  be  written  in  operational  form  or  in 
the  form  of  a Feinman  continual  integral* 


J 
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(7) 


where 

Do®  j]  *»(«)(  f ...  f j[  rfti(‘)exp  i i f 
{.-o  l J J {■  o 2 6 

In  the  first  of  the  equalities,  (7)  after  carrying  out  the  function  of  the 

operator  in  the  exponent  for  the  brace  which  follows  it,  one  must  place 

0,  A clear  presentation  of  the  solution  in  fora  (7),  in  several  cases 

is  suitable  for  investigation  (see  [12}  and  ), 


2,  Approximation  for  delta-correlated  fluctuations  of  a dieleotric 
oonstant  along  the  propagation  direction 

As  was  said  above,  field  u(x, fi)  only  funotionally  depends  on  pre- 
ceding values  of  f (£,/?).  However,  there  can  be  a statistical  link  between 
u(xt/fl)  and  subsequent  values  a^x',^'),  since  values  f(x',^')  when 
x' >x  are  correlated  with  values  of  c(£,  fi)  when  £<  x.  It  is  obvious 
chat  the  correlation  of  field  u(x,^?5  with  the  subsequent  values  of  I(x* ,/tf' ) 
is  marked  when  x*  - x ^ • where  1^  — is  the  longitudinal  radius  of 

correlation  of  £.  At  the  same  time,  the  characteristic  radius  of  correlation 
of  field  u(x,yfl)  in  a longitudinal  direction  has  a value  in  the  order  of 
x (see  for  example  [l]),  Therefore,  in  the  problem  examined,  there  is  a 
small  parameter  of  l^/x,  which  can  be  used  for  an  approximate  solution.  One 
can  place  l|a  0 in  the  first  approximation.  In  this  case,  values  of  fields 
u^i»  whnn  4 i<*  will  be  not  only  functionally  but  also  ntatistically 


independent  of  values  of  ¥(rjj.  fi  j)  when  > x,  that  is 

i i«o„  p«)h\.  tj)'y  - < ii  Pi))-  <(  n * *ty 

/ \ i ^ ' (8) 

(5,  < x,  ^ > x). 

By  using  the  feature  of  (8) , one  can  easily  find  equations  for  the  statis- 
tical moments  of  field  u(x,yj  ).  Let  us  indicate  this  on  an  example  <u>  . 
Let  us  use  equation  (4)  for  this.  By  averaging  it,  we  take  into  account  that 
in  the  seoond  term  in  the  right  part  the  value  e(rj,  fi)  in  the  exponent  is 
always  taken  when  values  of  ,i>J  , that  is,  statistically  independent  of 
the  seoond  factor  A pj  . Therefore,  when  averaging  (4),  these  factors 
can  be  averaged  independently)  x 


< «(*.  p)  > « « o(p)  ^cxp  | i ~ J p)  )>  + 

X k 

+ ^frf?^OXPf/2  P)  p)>’ 

0 I 


The  equation  obtained  la  closed,  since  it  does  not  contain  other  known  func- 
tions, apart  from  <u>.  In  a similar  manner,  one  oan  obtain  an  equation 

for  the  moment  of  the  arbitrary  ordert 

p„  p„;  pj O ~ < u(x,  pi) «(*,  p„)  X 

x u*(x,  p|), ....  «•(*,  p;„)>, 

For  this,  one  must  first  of  all  write  the  differential  equation  for  u(x,/#1), 
mi,  u*(x,  y0*m)  and  then  oonvert  it  into  an  integral-differential  oquation 
of  type  (4),  After  this,  by  using  (8),  one  can  carry  out  averaging.  An 
equation  for  M0t  n obtained  in  this  way  in  pi].  has  the  form 

|pf|,  !p/j)  “ M#,  »i((P/)i  iP/|) 


,r  n tn 

«/ ; p,)-^]t(h  p‘()  )V" 

o II  i i 

v x n tn 

2*  H j 1 \ exp  1*  2 j t!r'  ^ ^ p,)  ~ * ( * p<*  }y 
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where  (/>■*,  { is  the  aggregate  of  all  and  the  value  of 

Al2,«(|p/I,  (p/I)  - «u<Pl) «i>(pM)«u(pi) «o(pmj- 

la  designated  by  ^°n  n* 

3.  Equations  for  moment  of  a wave  field  In  a medium  with  a Gauss 
distribution  of  fluctuations  of  the  dieleotrio  constant. 

When  deriving  equations  (10),  only  the  oharaoteristlos  of  the  delta- 
correlation  of  fluctuations  of  the  dielectric  constant  along  the  direction 
of  wave  propagation  were  used;  the  distribution  law  itself  was  not  defined 
oonoretely.  Let  us  now  examine  a particular  case  of  the  Gauss  law  of  dis- 
tribution for  St  In  this  case,  the  statistical  characteristics  of  S are 
completely  described  by  the  correlation  function 

lh(x,  f;  x\  p')  ;'«(*,  f)i{x\  p')>. 

The  delta-correlation  condition  is  equivalent  to  substituting  Bg  by  an 
effective  correlation 1 

p;  p')  H*~* ')Mx,  P,  p'),  (11 

* 

here  A is  determined  from  the  equality  of  integrals  from  ^ and  B|  i 


raspeot  to  x't 


w 

A(x , p,  p')  f dx'Hx,  p;  p'). 


In  a case  of  statistical  uniformity  of  the  field  if,  which  we  shall  examine 
for  simplicity,  A * l\(p-  p ')»  that  is  A does  not  depend  on  x.  In  the 
oase  examined,  one  can  easily  find  the  value 


/esp 


and  after  calculating  it  reduce  equation  (10)  to  a differential  equation  [lij 
~ M„.  m ~ — { Aji  -f- ...  -}-  A„  — A,'  — ...  — AJJ  Mi,  m — 

v,  (12) 

~ ” Q(Pi.  •••!  Pm)  Mn,  m, 


Whare  S V,l(p, fl)  + 

/ i / i /-i  »--t 

m ;/i 

+ S (13; 

In  the  case  examined  of  Gauss  delta- correlation  fluctuations,  the  1 
aggregate  of  equations  (12)  and  also  their  interconnecting  equation  for  the 
characteristic  functional  of  a field  u can  be  obtained  by  another  method 
[10]  , vjhich  we  shall  illustrate  on  an  equation  for  <u>,  By  averaging  (2), 


we  obtain 


<«(*,  *)>  + **(«(.v,  t)u(x,  p)>  0. 


To  find  the  last  term,  we  use  the  formula  obtained  by  Purutsu  [29]  and 
Novikov  [3^ 

<-(D*[«i>  jdr'  (t(r)t(r');  <aA»|«] /»7(r')>,  (15) 

allowing  us  to  calculate  the  oorrelation  of  the  Gauss  random  function  t (r) 
«*>  a 0)  with  the  functional  R[i]  from  it.  Formula  (15)  can  be  proved, 
for  example,  by  expanding  R M into  a functional  Taylor  series, 

Slnop  the  solution  of  u(x,^)  of  equation  (2)  is  a functional  of  1, 
one  can  write 

<«(*,  p)c(*',  p’)>  />’.(*.  p,  p'), 

<*(*,  p )«(*,  p)>  \dx'  \ (I  p'  B,(x,  P,  p')\ 

S 1 

< \ (16) 

p /" 

By  substituting  instead  of  ^(x,^  x1  the  effective  correlation  func- 

tion (11),  we  obtain 


<~(*.  p)>  “ — frfp ‘A(p  - 

~ J '«*(*,  PV 

where,  by  integration  with  respect  to  x,  the  parity  of  the  ((-function  is 
taken  into  account,  as  a result  of  which  the  1/2  factor  appears.  By  sub- 
stituting the  averaged  ratio  (6),  we  obtain 

< i(.t,  f)  tt(x,  p)>  ■-  / >1(0)  < u(x,  f)  > , 

4 

and  equation  (14)  takes  on  the  fora 

(2tt^  + Al)<"(*,p)H'/V^0,<“(jr*p)>  ■"  °-  (17) 

Equation  (17)  coincides  with  equation  for  0,  obtained  from  (12), 

A complete  statistical  description  of  field  u(x,^)  can  be  obtained 
from  the  characteristic  functional 

UVJtvt'*]  - exp  {/  y d p («(.«,  f):'(p)  + «•"(.«,  p)<,!'<p))|  >• 

B(y  differentiating  fx  with  respect  to  x,  by  using  equation  (2)  and  ratio 
(6),  an  equation  for  )lf  x can  also  be  obtained.  It  has  the  form  \jo] 

0x  2k J r\  KYI  n v(p)  B r*(p) 

*r  A A <18‘ 

where  1 J 


;V/(p)  «.  ■£<(>)  - 

r<  i'(p) 


“<?) , ;; . 7* 

'*  i'  ip) 


Equation  (18),  as  can  easily  be  shown  flO,  15]  is  equivalent  to  the 
aggregate  of  equations  (12),  Equation  (18)  is  the  infinite  dimensional  analog 
of  the  Einstein-Fokker  equation,  and  in  this  respect  the  described  approxi- 
mation of  wave  propagation  in  a medium  with  Gauss  delta-correlation  fluctuations 
of  f can  be  called  diffused. 

Let  us  write  in  clear  form  the  equation  for  the  function  3 ^ ^ 

(\H(X’  R + 2 p j"'’(  *-  R -■  \ p)  > and  for  function  i4,A%  f|<  f,f  ^ ^ 

' u(\,  p,);/(.v,  p.i.''  '(a,  ?J)h" (a,  p’),  , resulting  from  (12),  when  n = m * 1 and 


n = m = 2: 


ac.  l /5s  b- 

~ - T - T WrfG*.  W9)  - ••'((')  ~ Ait), 

dx  k d A?  dp  4 

G=(0,  •?,  ?)  - u,(n  + ~ p)«j(*  - \ pj; 

--4  ...  L |Al A,  - a;  - a;i  id<Pi  - r,;>  + D(Pi  - ^ + 

•f  n<p,  - pi)  + d{?1  - ?;i  - D(h  - ?I)  - % ; - ?;)ig4, 

g.(o,  pi,  p...  p;,  p2)  - mpi)  woW  «i(p;>  «;<pi)- 


(19) 


(20) 


Equation  (19)  was  obtained  by  Dolin  [61  , and  Beran  [3ll  * Equation  (20) 
was  first  obtained  by  Shishov  [7j.  Equation  (12)  for  moments  of  arbitrary 
order,  apart  from  works  mentioned  above,  was  also  obtained  by  another  method 
in  a work  by  Chernov  [8]. 

As  Dolin  noted  [6],  equation  (19)  is  equivalent  to  the  so-oalled  narrow 
angle  approximation  of  a radiation  transfer  equation.  In  fact,  by  substi- 
tuting (19)  Gj(*f  ft.  })  , f R,  P), 

one  can  obtain  an  equation  for  function  for  J 

°.i  ■ .1  X/}j  ...  7 J + J*  - v.')./(.v.  K,  x'),  (21j 

where  7 ^ d x /(-/•)  — is  the  coefficient  of  extinction, 

- *- 

/(*)  “ -MM  — is  the  scattering  indicatrix,  'M*..''  . 7 1 — is  the 

three-dimensional  spectral  density  of  the  correlation  function  Be(r). 

Equation  (19)  for  Gg  can  be  solved  in  a general  case  for  the  arbitrary 
function  D { p)  and  arbitrary  initial  conditions.  As  for  an  equation  for  G;+, 
it  cannot  be  solved  analytically  and  one  must  resort  to  numerical  methods. 

Apart  from  equations  for  the  mean  values  for  the  product  of  fields  u(x,  /l>. 
u*(x, where  all  the  arbitrary  x are  identical.  Equations  can  also 
be  obtained  for  functions  <u(x., , , ...,  u*(x*  for  non-concurrent 
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1 

values  of  x [I4j,  The  boundary  conditions  for  these  equations  will  oontain 
functions  of  K for  concurrent  values  of  arbitrary  arguments. 

Let  us  also  note  that  presenting  the  solution  in  an  operational  form  or 
in  the  form  of  a oontinual  integral  (?)  when  using  suppositions  on  the  Gauss 
distribution  of  f and  the  delta-. correlation  along  an  arbitrary  coordinate 
makes  :\t  possible  to  calculate  < u(x,^)>  and  G2(x,R,^),  In  these  cases* 
it  is  also  possible  to  solve  equation  (12).  In  cases  when  the  solution  of 
equations  (12)  is  not  possible,  formulation  (7)  also  makes  it  impossible  to 
find  a clear  expression  for  Mn  m,  or  however,  it  i3  suitable  for  investi- 
gating the  asymptotic  behavior  of  moments  06J. 


L 


4.  A method  of  successive  approximations  and  a condition  for  the 
applicability  of  a diffused  approximation 


Let  us  now  oxamine  a more  general  method  for  obtaining  equations  for 
moments,  a particular  case  of  which  is  the  approximation  of  a diffused  random 
process.  Let  us  illustrate  this  method  by  an  example  of  an  equation  for  an 
•average  field. 

We  shall  reckon  that  the  fluctuations  of  f hive  a Gauss  distribution, 

df* 

but  the  correlation  function  De  (x,^v  x',^1)  is  not  changed  to  * . 

Let  us  average  equation  (2)<jand  to  find  < t we  use  the  Fui-utsu-Hovikov 
formula  (15).  As  a result,  we  obtain  the  expression  (16),  the  substitution  of 
which  into  the  equation  for  <'J ) gives 


( 2 tk?-  :•  A,  ) < « > + k'  \dx’  J d p7?(. r,  p,  p') 


" S«(x.  p)  \ 

" • p') 


(22) 


0. 


/u  //\ 

Equation  (22)  is  not  closed,  since  it  contains  a now  unknown  function  \ ~ N 


To  obtain  an  equation  for  this  function,  we  shall  act  on  equation  (2)  by  the 


1 2 


operator 


0 


m.- 


l«(*'  p')  when  x‘  «C  x and  w®  -hall  average.  For  finding 
value  of  y • (*.  p)  \ we  again  use  formula  (15)*  As  a result,  we 


obtain 


8*(*'.  p') 


I 2 Ik  --+■—)  + *>  (dx"  f tfp"B,(x,  p,  p")  > 

V ,«*  'Nl.^'.p')  / J 


y 


(33) 


y y_.  i_\  » 0. 

s »«(*\  p')Kv*.  P")  X 

Equation, (23)  is  again  dosed,  since  it  oontains  a second  variational  deri- 

• i- . 1 

vative  from  u.  Efcr  continuing  this  process,  one  can  obtain  an  infinitely 
linked  chain  of  equations.  A dosed  equation  in  approximation  of  a diffused 
random  process  is  obtained  when  converting  in  equation  (22)  function 
% to  Bjff  * <«x  - x')  A(^  -/>’),  since,  in  this  oase,  there  ia  a value 
<3m(x,  p )/3*(.v\  p')>  for  concurrent  values  x a x1 , which  according  to  (6) 


is  expressed  through  <u>.  Converting  B,  to  Bg  can  be  done  not  in 


the  first  equation  of  the  ahain  (22) , but  in  on?  of  the  subsequent  equations 

/ >»«(*,  p) 


For  example,  if  this  is  done  in  (22),  value 
under  the  symbol  of  the  integral.  It  can  be  obtained  from  equation  (6)  (in 


--  : -\  WjH  appear 

N««(*,P'W,P  y 


which  0X  is  converted  to  fl*) , if  one  acts  with  tne  operator  on  this  equality 
8 

— and  average 

y Vu(jc,  p)  v 


/_ 


\7(*.  f")w,  p'j  / 


Ht(x\  p')  / (24) 

By  substituting  this  expression  into  (23),  ws  obtain  a closed  equation  for 

/3u(*,  p)  \ 

/•  Aftor  this,  the  system  of  equations  (22),  (23)  can  be  solved, 
N3e(x',  p'r 

and  we  shall  find  a second  approximation  for  <u>. 

A second  approximation  is  much  more  accurate  than  an  approximation  of 
a diffused  random  process.  It  Is  quite  difficult  to  show  it  in  general  form, 
but  this  question  can  be  investigated  in  an  example  permitting  an  accurate 
solution.  As  an  example,  let  us  examine  equation 


1 


■J 

■ 
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(ill  “* 

2ikdx  + *1,(  v)',(Jc)  0 («(0)  * w0), 

which  differs  from  (2)  by  the  lack  of  term  u.  If  one  introduces  8,1*.  ■O 
” *?exp(-|*  — x'\ll),  % then 

i7(.V)  :v-  W0  f‘Xp(  • - |l(T  --  1 -|-  (?•"')], 

where  t •»  xjl,  n --  — /* <jJ,  i — is  the  radius  of  correlation  $ , 

4 1 

the  solution  of  the  first  and  second  approximations  respectively  have 

the  form  «i(x)  - «0exp(-»  p*), 

Hj(jt)  - (cxp(— pt)--  i«. exp( — t)l. 
t -p 

By  comparing  • u2»  and  u when  1,  one  can  easily  determine  that  the 
approximated  solutions  oan  only  have  a good  approximation  when  /*<  1 . Apart 
from  this,  funotion  u2,  as  opposed  to  u^,  at  aero  has  the  some  form  as  an 
accurate  solution.  Figure  1 shows  functions  of  u^,  ug  and  u when  ju  » C.2 
(•••  - u/u0,OOO  . u^/uq,  XXX-  u2/u0).  In  this  case,  the  difference  in 
the  Second  approximation  from  an  accurate  solution  does  not  exoeed  2.5$. 
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Equations  of  a second  approximation  can  be  solved  for  <u^  and  Gg. 

A comparison  of  the  first  and  second  approximations  makes  it  possible  to  find 
conditions  for  the  applicability  of  approximation  of  a diffused  random  prooess. 
Here,  it  appears  that  the  conditions  for  the  applicability  for  have  the 

form 

1,  *>tt  M > 1 , 

(25) 

where  -j  *»  — /(0)  — is  the  coefficient  of  extinction,  1 — is  the  radius 

•1 

of  correlation  of  a field  ¥ • At  the  same  time,  conditions  for  the  applioa. 
billty  of  a first  approximation  for  Oz(x,R,^)  have  the  form 

P\U\  1 f)  sC  M. 

where  — is  the  structural  function  of  a complex  phase,  calculated 

by  the  continuous  perturbation  method  ia 

It  must  be  stressed  that  conditions  (25)  and  (26)  are  practically  in- 
dependent, since  they  apply  limitations  to  various  parameters.  In  particular, 
it  can  be  shown  that  conditions  (2b)  are  fulfilled  wuen  condition  Vl^  1 is 
violated.  Let  us  state  that  condition  (26)  sets  limitations  only  on  the 
local  characteristics  of  fluctuations  of  If  and  therefore  can  also  be  written 
for  a turbulent  medium. 

Seditions  for  the  applicability  of  a parabolic  equation  (2)  are  inves- 
tigated in  work  £12]  , where  comparison  of  equations  solutions  (1)  and  (2), 
written  in  operational  form  (7),  were  done.  Here,  it  appeared  that  conditions 
for  the  applicability  of  a parabolic  equation  were  done  throughout  the  whole 
range  of  applicability  of  a diffused  approximation. 


5.  Solution  of  an  equation  for  the  function  of  the  mutual  coherence 
and  some  experimental  data 

Equation  (1?)  for  the  average  field  can  be  easily  solved.  Its  solution 
has  the  form 
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IP- 


t 


f 


< n(x,  f>)  > - «0(*,  p)o.xp|-  x j, 


k* 

where  uQ(x,^  ) — is  the  solution  where  there  are  no  fluctuations, T « - ,4(0)  — 
1»  the  ooefflolent  of  extinction. 

In  a general  oase,  one  can  also  solve  an  equation  for  function  Gg. 

This  solution  was  obtained  by  Brenner  j^2‘}  and  Dolin  [?3)  when  investigating 
the  equation  for  the  transfer  of  radiation  in  narrow  angular  approximations 
(21),  equivalent  to  equation  (19),  before  an  equation  was  obtained  for  Og. 
Later,  a similar  solution  was  investigated  in  works  ^13,  3^*  Whereas  in 
(19)  one  oarries  out  Fourier  conversion  on  variable  R,  which  is  not  included 
in  coefficients  of  the  equation,  we  obtained  a linear  differential  equation 
in  partial  derivatives  of  the  first  order,  which  can  U Midly  solved  bey* 
the  characteristic  method.  The  solution  has  the  following  form 

Gj(.v,  R , p)  =->  J dR'dp'  6».(0,  R •"  R\  p — ?')  X 


d \D\  p ~ P 


, t 


(2?) 


r ikR't'  ft*  f 

’■  ’ *! 

where  Qg(0,  R,^)  - is  the  initial  value  of  the  coherency  function  when 
x ■ 0.  Formula  (27)  can  also  be  used  for  partially  coherent  sources 
when  Qg(0,  R,^> ) do  not  have  the  fom 

Uo ( R -f  J 

For  a partial  oase  of  a plane  incident  wave  when  Q2(0,  R ,/))  • const  =* 
|uQ|2  , the  integral  along  R»  in  (27)  is  calculated  and  produces  rfOc/Wx), 
after  which  the  formula  takes  on  the  simple  form  oft 

G,(x,R.p)  |«.rvxp|  • Y |.  (2?a) 

If  one  examines  the  fluctuations  of  a dielectric  constant,  caused  by  turbulent 
fluctuations  of  temperature,  in  the  considerable  interval  of  wave  numbers  the 
three-dimensional  spectral  density  (H)  has  the  form 
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*!*.(*)  **  ' (*mni  *man)» 

where  A ■ 0.033  — is  the  numerical  oonstant,  cj  — la  tha  structural 
charaotaristio  of  fluctuations  of  a dialaotrio  constant,  depending  on  meteoro- 
logloal  conditions  [0*  In  this  oast,  function  D(^)  la  calculated  and 

#qU&1*  D(f)  w 2*  Jrf*  «»>.(*)[  I - ex  p(<«t)l  - A/CjpW. 

where  N * 36er*  aq<7 /6)/5  • 22^Q(1l/6)  ■ 1,46,  By  substituting  this  expres- 
sion Into  (27),  one  oan  calculate  both  as  a function  02 (x* t &)•  and  as 
average  intensity  of  a wavs  <I(x,  R)>  * 02(x,  R,  0). 


' ■ 0,9  -0,6 

. — ~r— w«»i  -hi 


«.i  Cf  'aS 

• i...  >.i.  «.  <u.  • ^ i . ..>.ui  J.a> 


au; 

*\>V 

"•Aft 

% 


Figure  2.  The  space  speotrum  of  the  coherency  function  of  the  seoond 
order.  Points  — experimental  data  according  to  [,3oJ.  the 
continuous  curve  was  plotted  on  the  basis  of  formula  U7) 
with  respect  to  the  effect  of  the  internal  scale  of  turbulence, 
described  by  parameter  -((cm  [)6]), 

Figure  2 shows  the  spaoe  spectrum  of  the  coherency  function  obtained  in 
work  [36],  On  the  same  graph  is  the  theoretical  dependence,  plotted  on  the 
basis  of  (27).  In  work  (j3?3  the  behavior  of  space-limited  beams  in  the 
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atmosphere  was  experimentally  studied.  In  Figure  3,  the  results  of  measure- 
ment of  average  intensity  are  oompared  with  a curve,  plotted  on  the  basis  of 
formula  (27), 


Figure  3 

Aa  for  the  equation  for  0^,  it  was  not  possible  to  solve  it  analytically. 
An  approximated  solution  of  equation  (20)  (with  respect  to  single  scattering 
in  the  sense  of  the  theory  of  radiation  transfer)  is  shown  in  [13,  1 5]  for 
turbulent  fluctuations  of  a dielectric  constant.  In  work  by  Dagkesamanskaya 
and  Shiahov  [j8]  there  are  results  of  a numerical  solution  of  this  equation 
for  the  correlation  funotion  of  a dieleotric  constant  in  the  form  of  a Qauss 
curve.  In  this  work  they  obtained  saturation  of  intensity  fluctuation,  which 
agreed  qualitively  with  results  described  in  [2  - 5)# 

Let  us  state  that  function  0^  also  describes  fluctuations  of  the  shift 
of  space-limited  beams  in  random  irregular  media  [13,  39]* 

6.  Amplitude-phase  fluctuations  of  a light  wave 

Let  us  now  examine  the  statistical  description  of  amplitude- phase  fluc- 
tuations of  a light  wave. 
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By  introducing  a oomplex  phase  of  wave  (jp  * =>  IS,  where  ^ * In  A/Aq  — 

is  the  amplitude  level,  and  3 — are  the  fluctuations  of  the  phase  of  a wave 
relative  to  the  phase  of  an  incident  wave  (x,  k),  one  can  write  the  equation 
for  the  continuous  perturbation  method! 

2/A^-t-  lj.<p  -j-  (Vj/r)'-’  4*  *’»  “ 0,  «»  exp(«p>. 

Aoourate  solutions  of  equations  (2)  and  (28)  are  equivalent.  In  the  first 
approximation  of  the  continuous  perturbation  method  in  equation  (28),  the 
term  (V.;.?)11  is  omitted,  that  is,  the  equation  examined  is 


2/A  2ft.  + Aj.<pfi  + A»«(*,p) -•<>. 
ox 


In  this  case,  the  solution  is  a Qauss  random  field  and  the  statistical  pro- 
perties of  amplitude  functions  are  described  by  the  parameter  a? »« < ^/a*,  P)>  [ll. 
Here,  the  validity  condition  of  the  first  appro ximatinn  of  the  continuous  per- 

p 

turbatlon  method  for  amplitude  funotion  is  the  condition  a 1,  In  the 

p 

range  «r  Q ^ 1 (which  is  oalled  the  range  of  strong  fluctuations)  one  must 
study  the  full  equation  (28). 

Let  us  state  that  the  diffused  approximation  for  equation  (2)  dues  not 
impose  limitations  on  ths  amplitude  fluctuations  and,  consequently,  equations 


for  momsnts  of  fitld  u(xl/(d)  (12)  ars  also  valid  in  the  range  of  strong  fluc- 
tuations of  amplitude. 

Statistical  properties  of  the  solution  (28)  can  also  be  described  in 
diffused  approximation.  However,  on  the  strength  of  the  nonlinearity  of  the 
equation  (28),  equations  for  moments  of  function  are  incomplete. 

By  separating  the  imaginary  part  in  (28),  we  obtain  equations  for  the 
level  of  amplitude  'JC  (x,^>)  and  wave  intensity  I = exp(2^)t 

2/e --  + Aj S + 2v,/v  ■ S - • 0;  (30) 

ox  x - 
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I 


f 


r 

f 


(31) 


Experimental  study  of  the  distribution  of  prob&bilitiss  for  amplitude 
level  [2]  showed  that  this  distribution,  both  in  the  weak  range  and  In  the 
range  of  strong  fluctuations  is  olose  to  a normal  distribution,  although 
"tails"  of  this  distribution,  at  the  present  time,  have  been  poorly  studied 
and  they  also  determine  the  distribution  of  probabilities  for  wave  intensity. 
As  for  fluctuations  of  the  angle  of  arrival  of  a wave  at  the  observation 
point,  linked  with  fluctuations  of  ^ S,  they  are  well  described  by  the 
first  approximation  of  the  continuous  perturbation  method  . 

For  a plane  lnoldent  wave,  boundary  conditions  for  (30),  (31)  are  con- 
ditions /A  f)  0,  /((),  p)i- 1 1 and  the  solution  of  these  equations  will  be 
similar  to  random  fields  in  plane  x « oonst. 

By  averaging  equation  ( 31 ) , we  obtain  the  ratio 

(32) 


expressing  the  law  of  conservation  of  energy  for  the  examined  problem, 

Gy  multiplying  (31)  by  ^ and  averaging,  we  obtain,  by  allowing  for  (32) 
and  apaoe  homogeneity,  the  ratio 

< P>/(*.  p)  > - “ J fM  < rtVjSU.  9)  y ( 33) 

0 

Values  in  the  left  part  (32),  (33)«  are  determined  by  the  single-point 
distribution  of  probabilities  for  field  and  the  value  in  the  right 

part  of  (33)  is  linked  with  the  correlation  of  value  and  ^S, 

If  the  single-point  distribution  of  probabilities  .for  % is  of  the  Gauss 
typs,  then  aooording  to  ("12)  it  will  be  the  expression 

<■/(*.*)>'■-  (34) 

2 

where  o — ia  the  dispersion  of  the  amplitude  level,  in  this  case  <//>  -»*. 
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If,  whan  calculating  the  anplltuda  fluctuations, tha  supposition  on  tha 
possibility  for  using  tha  axprassion  from  tha  first  approximation  of  tha  con- 
tinuous perturbation  method  for  Yi  S is  valid,  tha  right  part  of  (33)  oan 
be  calculated  In  a diffused  approximation  tyl]  , and  wa  obtain 

**  In  °o  (=?»  1).  (35) 

In  this  way,  tha  result  of  tha  oombinad  use  of  suppositions  on  tha 
normality  of  a single-point  distribution  of  probabilities  for  amplitude  level 
and  tha  possibility  of  substituting  Vi  S into  tha  expression  from  tha 
first  approximation  of  tha  continuous  perturbation  method  is  expression  ( 35)  • 
whioh,  however,  does  not  ooinaide  with  experimental  data  relative  to  tha 
behavior  of  depending  on  cr*Q  ^ " **.  Fuller  information  on  the  ampli- 
tude-phase fluctuations  of  a plane  light  wave  oan  be  obtained,  whioh  is  limited 
by  the  validity  of  approximation  of  gaometrioal  optics* 


II.  Approximation  of  the  Geometrical  Optics 

1.  Diffusion  of  Beams  in  a Medium  with  Random  Irregularities 
From  the  striotly  formal  point  of  view,  the  following  is  reduced  to  tht 
Einsteirv-Fokker  equation.  Let,  j ; ^(s), I„(s))  satisfy  the  system  of  dy- 
namic equations 


(t  Ms) 
ds 


«*  *)  + /((?.  s), 


(36) 


where  v^,  s)  — are  the  determinate  functions,  and  f^,  o)  — are  the 
random  functions  (n  + 1)  of  the  variant  having  properties 

a)  f^(£,  s)  — a Gauss  random  field  in  (n  + 1)  -dimensional  space  (£,  s), 

b)  </,<?.  s)>  0, 


o)  Oitvi  ')/*.(?  I *’) 


2?(*  - *')F,Al  V,  *). 


(37) 


In  this  oase,  the  probability  density  for  solving  system  (36),  that  is, 


function 


?Ax)  -<«(!(•)  -.v)> 


(hep*  ^ (■)  — is  solution  of  (36),  corresponding  to  the  specific  realisation 
of  f(<£,  s),  and  averaging  is  done  on  the  population  of  all  realisations  of 
f)  satisfies  the  Einstein-Fokker  equation! 

^+^IK(^S)-M,(*,S)|P.W|_  (j8) 

— dx~dx~  ^X'  *' 

Here,  the  notation  was  introduced  * 

OX i y s4  ,v 

and  summation  is  don*  on  the  reourrent  indices. 

These  facts  are  well  known  (see,  for  example,  [42,  43]),  For  praotioal 
use,  a suitable  deduction  for  equation  (38),  based  on  the  use  of  the  Furutsu. 
Novikov  formula  (15)*  is  shown  in  work  [2*3  . In  the  same  work  there  is  a 
method  for  finding  corrections  to  equation  (38),  allowing  for  a finiteness 
of  a radius  of  correlation  of  field  f along  the  0 axis,  which,  by  the 
same  token,  allows  one  to  determine  limits  for  the  applicability  of  the  equa- 
tion (38)  for  aotual  physioal  systems. 

Let  us  turn  to  equations  of  beams.  Normally,  beam  length  1 is  used 
as  an  independent  variable  s.  By  accepting  1 as  the  independent  variable, 
we  shall  write  equations  for  beams  in  a form 

where  * In  n,  n — is  the  refractive  index.  If  one  introduces  6-veators 
g « jr,  t|,  v «»  (t,  OJ,  /■»  |0,  a],  where  l7i  “ (5<*  “ V*) . , the  system  of 

equations  (39)  can  be  written  in  form  (36),  Conditions  (3?a,  b)  can  be 
accepted  without  reservations.  However,  at  for  the  transfer  to  approximations 
of  the  correlation  function  for  f by  using  the  <f-function,  we  meet  insur- 
mountable difficulties.  The  fact  is  that  •:)  <<,(;)  and  does  not 

depend  on  s * 1,  Fbrmally,  one  can  reckon  that  the  function  g\  s,  s') 
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(in  this  oass,  not  depending  on  s,  s')  has  according  to  variable  s - s'  an 
infinite  interval  of  correlation  for  with  as  large  an  increase  as  possible  in 
the  value  of  s - s',  but  with  ^ ^ ■ fixed,  this  funotion  does  not  diminish. 

In  this  way,  it  is  not  possible  to  write  the  EEF  corresponding  to  the  dynamic 
system  (39)*  However,  one  can  show  equations  of  beams,  by  taking  the  indepen- 
dent variable  to  coordinate  x.  If  one  looks  for  the  equation  of  a beam  in 
form  Ra  ■ Rx(x),  where  Rj_  *|y,  ij  — is  the  transverse  shift,  instead  of  (39) ( 
we  shall  have  a dynamic  system 

_7jl 


5**-/  r 


d 2 t ’ ?} 

dx 


VT-*l  <«>) 

where  t j,v>  ,f](  ,,,  (fia,  • However,  we  notice  that  in  this  form  the 

equation  of  beams  aan  be  used  only  until  the  first  point  of  turn,  where  the 
denominator  I 1--^.  is  reduced  to  aero.  It  follows  that  in  a statistical 
problem  (40)  can  only  be  used  in  a range  where  the  probability  of  negative  fx 
is  small,  that  is,  in  the  range  of  narrow  angular  deviations  of  beam.  Since 
in  this  range  1 1 — /*•!.  , instead  of  (40)  we  obtain  an  approximate  system 
of  equations  for  beams  in  narrow-angle  approximations j 


dx  Rl('X^  **  • x)‘ 


dx 


(41) 


For  system  (41)  s - x is  already  included  in  the  number  of  arguments  f, 
and  therefore,  here  one  can  transfer  to  the  EEF,  The  corresponding  EEF, 
according  to  ( 38)  has  the  form 


dx  > L ok' 


D 


0V>X 


(42) 


where  D — is  the  coefficient  of  diffusion,  whioh  ocours  when  calculating 
Fkl  and  in  a model  of  statistically  uniform  and  isotropic  fluctuations  of  ^ » 

W 

D ■ rs  \ !•(*). 


(43) 
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Hera,  <piH)  — is  the  three-dimensional  spectral  density  of  the  correlation 
function  ^ . Expression  (43)  for  D ooincides  with  those  shown  in  works  fid  . 
213.  Equation  (42)  is  easily  solved,  and  its  solution,  corresponding  to  the 

t 

initial  condition  'A fix)  has  the  form  of  a Qauss  distri- 


(44) 


bution  with  moments  < R±l  {,}  R , i(x)  s A D ^ 

<Rxt{x)  ".;.*(*)>  - 

*<*)  > ■■■■  2/nM,A*. 

These  expressions  are  well  known.  On  the  basis  of  system  (41)  it  is  also 
easy  to  obtain  the  longitudinal  correlation  function  of  beams  displacement 
C?4]  j 

_jg,i LM.&R  ,,  ( , + 1 At  + « fsL. 

V'<^(.«)><7J\l.v')>  \ 2 / Xml,.  (45) 

Let  us  now  examine  the  problem  of  simultaneous  diffusion  of  two  beams.  In 

this  case  we  have  the  following  dynamic  system  of  eighth  orderi 

(I  R 1 t _ if , a 'i  | A'  1 ,,  ,u 

A/A  * 4‘  "7.v“  ' <>/{.,  ' 

i = 1,  2 — number  of  beams.  The  corresponding  EEF  for  the  system  (46), 
acoording  to  ( 39) , has  the  form 


d n in  n - a - "A  1 - d/)‘- 

“i'j  I ‘III  -ill)  1 * 1 I -----[ 

(At  ('/f,  I <>/v’  } 


(4?) 


dJl\ 


o . 


0 • , ] 


ir  (A’,.  • A\,)/Jc, 


where 


1 rtvlfp) 

2 J.v . 0\ 


U”, • ; — .v,  y,  « r 
<4(p)  2-  I ff  xp).  p - A’  ; | !i  ■ 


In  equation  (47)  one  can  introduce  new  variables  p ■ v ft... ; A'ii), 

/ — "til  1 v,  “ i l " a) , After  transferring  to  new  variables,  it  is 

possible  to  carry  out  integration  on  Ft  and  C and  obtain  an  equation  for 
the  function! 

/)  < ;'( A*  | — p)  , /) 
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describing  the  relative  diffusion  of  two  beans, 

Of  Ol-Jll. 


ox 


Here 


IM?)  21  D3.J—  r^(p)|  - 2*  f d**,*  .«!•(«)  (I  -cos*?]. 


(48) 


(49') 


If  1Q  — is  the  radius  of  correlation  of  gradients  of  the  refractive 
index,  when  p>/0  D ,,,(?)  2/>3v.  This  ratio  shows  that  then  there  are  large 

initial  distances  between  beams,  in  comparison  with  1q,  their  relative 
diffusion  occurs  with  a two-fold  diffusion  coefficient,  which  corresponds 
to  the  independent  diffusion  of  each  beam.  The  joint  distribution  of  pro- 
babilities in  this  case  can  be  reckoned  as  Gauss. 

In  a general  case  it  was  not  possible  to  solve  equation  (48) , It  is 
only  obvious  tiiat  when  there  is  a variable  diffusion  coefficient  (49)  the 
solution  is  not  a Gauss  distribution.  This  can  indicate  [24  J that  the  mean 
square  of  the  distance  between  beams  increases  proportionally  to  the  cube  of 
the  distance  traveled,  but  the  coefficient  in  this  formula  is  different  in 
small  and  reciprocal  distances.  These  two  ranges  are  separated  from  each 
other  by  a small  range  of  exponential  increase,  the  beginning  of  which  coin- 
cides with  the  beginning  of  the  range  of  strong  fluctuations  of  wave  intensity. 

Let  us  examine  the  limits  for  the  applicability  of  a diffused  approxima- 
tion for  beams  shown  above.  As  was  already  said,  the  EEF  for  diffusion  oi‘ 
bsamB  could  only  be  substantiated  in  narrow-angle  approximation.  Hence,  according 
to  (44)  the  condition  arises 


or 


/;• 


x 


(' 


(50) 


where  — is  the  dispersion  of  fluctuations  of  the  refractive  index.  This 
condition  imposes  a 3mal 1 limitation  on  the  transverse  displacement  of  the 
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beams  x-.  An  examination  of  corrections  linked  with  the  finitflness  of 

the  longitudinal  radius  of  oorrelatio n,  does  not  lead  to  new  limitations. 

Let  os  state  that  in  another  extreme  case  (when  the  beam  goes  far  off 
into  a heterogeneous  medium)  there  is  an  isotropic  distribution  of  probabilities 
for  the  beam  in  space  similar  to  the  Boltzmann  distribution  in  statistical 
physios  0*4] , 

2.  Amplitude™phase  fluctuations  of  a light  wave  in  random  heterogeneous 
media 

Let  us  now  examine  the  amplitude-phase  fluctuations  of  a plane  light 

wave. 

In  an  approximation  of  the  geometrical  optics,  the  system  of  equations 
for  the  amplitude  level  and  wave  phase  takes  on  the  form 

+ 7 V ' a>.7.  — A S(x,  t)  (,:((),  p) « O), 

. + * 7frp)  (.5(0.  P)  »•  0).  ' C51) 

dx  ‘2k  2 

The  astern  of  equations  (51)  was  investigated  in  works  [45  - 4?,  2^, 

Solution  of  an  equation  for  *)C  can  be  written  in  the  form  [4,5p. 

X 

P)  - 5:p0,  (52) 

where  R(x’)  R (x,  x’\  p are  the  transverse  displacements  of  beams,  determining 
the  beam  tube,  which  satisfy  the  equation 


/?(*,  f)  - p - j J <n\  S(;,  R(x,  •;  p)). 


Here,  on  the  strength  of  the  .supposed  singularity  of  beams, 

R(h,  U /?(*.  «r,  f ))  «(x,  S;  P).  (5*i 

The  paths  of  beams  arriving  at  the  observation  point  (x, p)  are  determined 
by  the  value  Vj.  3 the  equation  for  which,  in  its  turn,  has  the  form 


Ratios  (52)  - (54)  aro  also  valid  when  there  are  diffraction  effects. 
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+ -£(vaSv*  )v.LS  =--  -*•  V t *(*!  p)-  (55) 

"he  eolution  of  equation  (55)  can  also  be  written  in  the  form  of  an  integral 


with  respect  to  the  same  bearni 


X 

VxS( x,  f)  - J j*  d 8(S,  /?( X,  5;  f)). 


By  substituting  (56)  into  (53).  one  can  disoard  the  wave  phase  from  examina- 
tion, and  we  shall  arrive  at  an  equation  for  the  path  [25] 


•r*  »(*',  k(x’)) 


with  limiting  conditions 


ot  , dli(x')  i n 
#(x)  p.  ~ v;H  “0i 

Here,  the  expression  (52)  for  the  amplitude  level  assumes  the  form 


/(*.  p)*> 


1 r d%  [<h, r(x.  v.  f))  x 
■1  j .1  «p/>p* 

0 « 


dty  f>  Aft  v,  «; 


In  this  way,  both  the  amplitude  level  and  the  phase  gradient,  determining 
the  angle  of  arrival  of  the  wave  at  the  observation  point  for  a specific 
realisation  of  field  l(x,  p) , is  determined  by  the  single  "dynamic"  equation 
(57)  with  limiting  conditions  (58) . Beams,  determined  by  equation  (57) 1 vary 
from  beams  in  Section  1 since  the  limiting  conditions  for  them  are  set  in 
different  space  points  (when  x*  * 0 and  x‘  « x).  Hence,  the  propagation 
of  beams  described  by  aquation  (57)  with  conditions  (58),  does  not  have  the 
nature  of  a diffused  random  process  (even  during  the  condition  of  6 -correlation 
of  field  ?(x,/>)  along  the  x axis),  and  the  probability  density  for  these 
beams  is  the  arbitrary  probability  density  for  a whole  aggregate  of  beams  in- 
vestigated in  Section  1 of  this  chapter. 

Beams,  determined  by  equation  (57)  with  conditions  (58),  describe  fluc- 
tuations of  the  beam  tube,  which  determine  fluctuations  of  wave  intensity  at 
the  observation  point.  In  work  0*8]  there  is  a calculation  of  fluctuations  of 
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the  sectional  area  of  a bean  tube  when  one  can  uso  a successive  approximation 
method  for  equation  (57) • 

Let  us  note  that,  as  was  shown  in  work  ^7] , both  the  solution  itself 
of  the  equation  for  •%,  and  all  the  statistical  characteristics  of  it  are 
determined  by  solving  the  aquation  for  the  wave  phase  S(x,  p)  ssifij, 

J,  The  applicability  of  diffused  approximation  for  describing  ampli- 
tude-phase wave  fluctuations 

An  equation  for  the  amplitude  level  can  be  written  in  the  form 

(/<0,  p)«l),  (60) 

where  /*»  exp(2/)  * itix,  pUt'^x,  p>/|  i/0  )*>  ““  ia  the  w*ve  intensity.  On  the 

etrehgth  of  the  space  homogeneity  of  field  l?(r)  random  fields  u(x, p).  I(x, p), 
S(x,p)  will  be  homogeneous1  fields  in  plane  x * const.  By  using  equations 
(51)*  (60)  one  can  obtain  an  expression  for  the  single-point  intensity  corre- 


lation and  wave  phase* 


on  the  other  hand,  the  rutin  takes  place 

Pu  Pi)  'Pi  t ■:  < /((*•  .>.)•  (vrS')’|>,  (62 

where  Gu.  p,.  p.)  ■■■■■  <//(.v,  p,)wix,  p.)  > is  the  mutual  coherency  function 
examined  in  the  first  ahapter.  In  an  approximation  of  geometrical  optica, 
the  expression  (62)  is  simplified  and  takes  on  the  form 

ViVtG(*>  pt,  P«)j^f-  </(Vj S)*>.  (6< 

The  left  portion  (63)  oan  be  calculated  in  an  approximation  of  the  diffused 
random  process*  In  this  approximation,  function  0 for  the  plane  wav*  is 
described  by  formula  (2?  a),  and  < *(.v,  p> /(.v.  p)>«0.  • Consequently,  in 

an  approximation  of  geometrical  optics  and  diffused  approximation^ 


Let  us  state  that  formulae  (64)  are  also  valid  in  the  cass  of  non-Qauss 
fluctuations  of  t , 
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(64) 


</$>«-£aa/)(0)*’, 

16 


» 1 

where  d**'®«(x)*  ■»  — < (vxS«)*  > — is  ths  dispersion  of  the  angle  of 

arrival  of  a wave  at  the  observation  point  in  the  first  approximation  of  a 

continuous  perturbation  method  (the  parameter!  characterising  the  intensity 

, * 

of  phase  fluctuations  6^] .,  S6(x,  p)  - ~ J p)), 


The  correlation  </$>«■  L < «w*|nw/w*  ^ is  determined  by  the  single- 
point distribution  of  probabilities  of  field  u(x,  ^0)  and  generally  speaking, 
is  conditioned  by  all  momenta  of  u(x,^).  The  expression  (64)  for  this 
correlation  was  calculated  in  a diffused  approximation. 

By  following  Section  4 of  the  first  chapter,  one  can  obtain  conditions 
for  the  applicability  for  formulae  (64),  whioh,  to  some  degree,  will  charac- 
terise the  condition  c f applicability  of  the  diffused  approximation  for  a 
single-point  distribution  of  probabilities  in  an  approximation  of  geometrical 
optics.  These  conditions  for  a turbulent  medium  will  have  the  form 


»**»!>  T«l.  (65) 

where  the  value  of  ^ characterises  the  internal  scale  of  turbulence. 

1'n  this  way,  when  carrying  out  conditions  (65)  calculations  for  the  am- 
plitude-phase .fluctuations  In  an  approximation  of  the  diffused  random  process 
in  a geometrical-optics  approximation  are  valid.  Conditions  (65)  coincide  with 
conditions  (50),  describing  ths  applicability  of  a diffused  approximation  for 
beams  in  a turbulent  modium,  which,  as  was  shown  in  the  previous  section,  also 
determine  all  the  statistical  characteristics  of  a light  wave. 

Let  us  try  to  analyr.e  the  results  obtained  from  the  point  of  view  of  their 
applicability  for  the  propagation  of  light  in  a turbulent  atmosphere. 
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A parabolic  equation  (qua si-optical  approximation)  was  suggested  as  a 
basis  for  examination*  This  equation  is  only  valid  when  one  can  disregard 
bade  scattering*  Pbr  light,  scattered  in  turbulent  fluctuations  of  the 
refraction  index,  angles  in  which  the  basic  part  of  the  scattering 

e k 

field  is  included,  have  an  order  of  10  v - 10  * Since  fluctuations  of  the 

dielectric  oonstant  in  the  atmosphere  are  also  very  small  (St — 10~CJ  , one 

oan  almost  always  disregard  back  scattering  of  light  in  turbulent  hetero- 
geneities* This  condition  is  only  expected  to  be  violated  in  the  millimeter 
wave  band,  and  even  in  this  case  deviations  are  very  small. 

On  the  other  hand,  the  small  value  of  the  scattering  angle  $ causes  a 
strong  interaction  of  scattered  waves  among  themselves  and  with  the  incident 
wave.  This  effect,  in  apite  of  the  small  value  of  & l also  causes  strong 
field  fluctuations.  As  a result,  as  was  said  above,  calculation  methods, 
based  on  one  or  another  form  of  the  perturbation  method,  are  invalid.  The 
suggested  theory  uses  the  small  value  of  another  parameter  — the  characteristic 
longitudinal  scale  of  heterogeneities  compared  with  all  other  longitudinal 
scales  of  the  problem.  Here,  there  is  an  analogy  with  the  nonequilibrium 
kinetic  theory  of  gases,  where  the  small  value  of  the  interaction  time  of 
molecules  (in  our  case,  the  longitudinal  coordinate  plays  a time  role)  is  used 
for  separation.  As  a result,  in  spite  of  the  fact  that  conditions  for  the 
applicability  of  the  theory  occur  and  contain  limitations  on  the  sizo  of 
fluctuations  of  St  * these  limitations  are  so  small  that  they  are  carried  out 
in  an  actual  atmosphere  even  for  relatively  long  routes.  Unfortunately,  due 
to  the  expansion  of  the  applicability  range,  one  must  oontend  with  a relatively 
complicated  theory  and  the  complication  of  its  results,  that  is,  one  must 
have  recourse  to  numerical  methods  of  calculation  (for  example,  for  function 
3^)*  Apart  from  this,  since  the  proposed  theory  operates  with  coherency 
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functions  of  different  orders,  complications  occur  when  phase  fluctuations 
must  be  calculated.  However,  the  basic  difficulty  with  which  one  must  con- 
tend in  the  continuous  perturbation  method  — of  not  being  able  to  describe 
a light  field  in  the  strong  fluctuation  range  — is  eliminated  in  the  new 
theory.  Apparently,  in  all  oases  when  the  statistical  characteristics  of  a 
light  field  in  a turbulent  atmosphere  which  interested  us  could  be  expressed 
by  coherency  functions,  the  suggested  theory  almost  always  gave  us  an  accept- 
able result. 
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